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Abstract 

The sum-over-states (SOS) polarizabilities are calculated within approximate mean-field electron 
theories such as the Hartree-Fock approximation and density functional models using the eigen- 
values and orbitals obtained from the self-consistent solution of the single-particle equations. The 
SOS polarizabilities are then compared with those calculated using the finite-field (FF) method. 
Three widely used mean-field models are used: (1) the Hartree-Fock (HF) method, (2) the three 
parameter hybrid generalized gradient approximation (GGA) (B3LYP), and (3) the parameter- 
free generalized gradient approximation due to Perdew-Burke-Ernzerhof (PBE). The comparison 
is carried out for polarizabilities of 142 molecules calculated using the 6-311+- |-G(d,p) basis set at 
the geometries optimized at the B3LYP/6-311G** level. The results show that the SOS method 
almost always overestimates the FF polarizabilities in the PBE and B3LYP models. This trend is 
reversed in the HF method. A few exceptions to these trends are found. The mean absolute errors 
(MAE) in the screened (FF) and unscreened (SOS) polarizability are 0.78 A 3 , 1.87 A 3 , and 3.44 
A 3 for the HF, B3LYP and PBE-GGA methods, respectively. Finally, a simple scheme is devised 
to obtain FF quality polarizability from the SOS polarizability. 

PACS numbers: 31.15.Ar, 31.70.Hq, 34.70. +e, 84.60.Jt, 87.15. Mi 
Keywords: polarizability, sum over states, finite field 
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In this paper, we examine sum-over-states (SOS) and finite-field (FF) polarizabilities for 



Hartree-Fock (HF), a hybrid density-functional B3LYP 1, 2], and a generalized gradient ap- 
proximation (GGA) pure density-functional, Perdew-Burke-Ernzerhopf (PBE)|3|], for a large 
set of molecules. These three are popular independent-particle models that are routinely 
used to study the electronic structure of molecules. We discuss the relationship between the 
SOS polarizabilities and the one-electron eigenvalues of these mean-field methods. We then 
devise a simple empirical scheme to obtain finite-field quality estimate of the polarizability 
from the HF and PBE models. We also investigate the differences in the predictions of 
polarizabilities from SOS and FF methods in the above popular single-particle methods. 

When a molecule is placed in a static electric field F, the electronic charge of the molecule 
redistributes. This response is characterized by the polarizability of the molecule. The 
molecular energy E in the presence of a small external field can be expressed as a perturbative 
series in power of F as 
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where, Eq is the total energy of the molecule in the absence of the external electric field, 
fi is the induced dipole moment, a is the linear polarizability, and (3 is the first hyper- 
polarizability; the indices i,j,k denote the Cartesean components. The calculation of the 
linear polarizability a^, which is the topic of the present manuscript, requires calculation 
of the second derivative of the total molecular energy with respect to applied field. The 
ctij terms can be calculated using perturbation theory or by finite difference technique by 
using the total energies calculated at various field strengths. A summary of various methods 
for practical calculations of polarizabilities can be found in several review articles^, [5]. 
Using standard time-independent perturbation theory, one can obtain the well-known SOS 
expression (see, for example Bartlett and Sekino in Refo). The SOS expression for a diagonal 
element of the polarizability tensor is given by 



Here, Eq is the ground-state electronic energy and Ek is the energy of the k-th excited state 
and and \l/o is the ground-state wavefunction and is the wavefunction of the k-th excited 
state. The matrix element (^ol^l^fc) is the z component of the transition dipole moment. 
Thus, knowledge of the excitation energies and the transition dipole moments can be used 
to compute the polarizability using the SOS method. The summation in the SOS expression 
(Eq. (131) is over all excited states. One major advantage of the SOS method is that the 
transitions that make significant contributions to the polarizability can be identified. If 
the perturbation is small then one expects the largest term in this expression to be the first 
because it has the smallest energy denominator. In mean-field theories this energy difference 
is often approximated by the difference between the energy of highest-occupied molecular 
orbital (HOMO) and the energy of the lowest unoccupied molecular orbital (LUMO). In the 
following we call the exact first energy difference, E% — Eq, and approximations to it the 
HOMO-LUMO gap. 

Exact polarizibilities can also be obtained using small finite fields. In the finite field (FF) 
method the molecular Schrodinger equation is solved for different directions and magnitudes 
of the applied electric field and finite differences are used to calculate the components, o^, 
of the polarizability tensor. The polarizability matrix elements can be obtained using a least 
square fit to the total molecular energies or by using a suitable numerical approximation to 
obtain the second derivative. 

The methods for the calculation of the molecular energy by solving Schrodinger equation 
can be broadly sorted in two classes: the traditional quantum mechanical methods such as 
the HF method and beyond, and and models based on density functional theory (DFT). In 
this article, we perform a comparative study of polarizabilities computed using the SOS and 
finite-field methods. The comparison is carried out for polarizabilities calculated within the 
HF approximation and two DFT models. 

Within mean-field theories such as the HF or DFT, the SOS equation [3] can be written 

as 

< ^m\Pi\^n >< 4>n\Pj\?Pm > 
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Here, the {ip} and {e} are the self-consistent molecular orbitals and eigenvalues, respectively. 
The summation in this expression runs over all occupied and virtual states. In practice, 
however, only a finite number of virtual states are included since the contributions from 
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very high energy states are negligible. Hereafter, we shall refer to the Eq. (j3J) as the sum- 
over-states SOS expression (and not Eq. (j3J), which is over fully correlated electronic states 
) and the polarizability calculated using this equation as the asos ■ 

The single-particle eigenvalues in Eq. (J3J) can be viewed as approximations to the excita- 
tion energies of the computational model. This expression treats each orbital as responding 
independently to the applied electric field, but as the orbitals respond to the applied field 



the self-consistent or mean-field component of the Fock and KS potentials also changes[6| 
Eq. (J4]) uses only the unperturbed eigenfunctions ipm and eigenvalues. The additional field- 
induced electron redistribution, or screening, is missing in the computation of polarizabilities 
by the SOS method. The asos are called the unscreened polarizability. 

The self-consistent treatment of field-induced polarization that includes the screening ef- 
fects can be computed by the FF method. In this method a term (-F ■ r) that represents 
interaction between the electrons and the applied field is added to the zero-field, one-body 
Hamiltonian, and the single-particle equations (Hartree-Fock or Kohn-Sham) are solved 
self-consistently. The self-consistent procedure takes into account the screening effects or 
field-induced polarizations. The polarizability is then calculated by numerical differentia- 
tion of the energy or dipole moment obtained at various field strengths in limit of vanishing 
field jl]. As a result of the complete description of the first-order screening effects, the cuff 
polarizabilities generally are found to be in good agreement with experimental values. Equiv- 
alently, the energy derivative required for the polarizability can be computed analytically 
by differentiation of the field-dependent Hamiltonian 0, ^]. The latter approach is called 
coupled HF method. 

The lack of an appropriate description of the screening effects in the asos calculations 
generally results in asos overestimating the experimental polarizabilities and op^in density- 
functional methods. Early calculations of Stott and Zaremba showed that the inclusion 
of screening by self-consistent calculation of polarizabilities of atoms reduces its value by 
roughly 40% [10]. Later, Mahan also noted that a self-consistent treatment reduced the linear 
polarizability by about 40% but the next two higher (hyper) polarizabilities are reduced by 
2% only 6|. Similar conclusions were drawn by Pederson and Quong, who in their calculation 
on the C 60 fullerene found asos to be roughly three times larger than the screened apF, 
with the latter value being in good agreement with experimental measurements Jll ] . These 
works show that the inclusion of screening effects is necessary for a satisfactory estimation of 
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polarizability. The appropriate treatment of screening effects using finite-field or equivalent 
methods requires full self-consistent solution cannot take advantage of the full unperturbed 
point-group symmetry. This lack of symmetry can make polarizability calculations computa- 
tionally intensive for large systems such as carbon fullerenes, quantum dots, finite nanotubes 

e schemes to correct asos f° r the screening effects have been 
Q. Recently, Gueorguiev, Pacheco and Tomanek have used a 
similar scheme to calculate polarizabilities of large carbon fullerenes within the tight-binding 
method. [15]. 



etc. For this reason, sim 
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The orbital eigenvalues that are required in the SOS method (Eq. (Tj0)) have different 
interpretations in HF theory and DFT [l^]. In HF theory, Koopmans theorem states that 



the occupied orbital eigenvalues are the negative of unrelaxed ionization energies. 



18|. On 



the other hand, in DFT only the highest occupied eigenvalue has a physical interpretation 



as the negative of the first ionization potential [19]. The DFT eigenvalues are the derivatives 
of the total energy with respect to orbital occupancies The practical DFT applications 
require approximations to the exchange-correlation functionals that model the exchange- 
correlations effects that are quantal in origin. Due to the nature of most of these approximate 
functional forms, the effective potential that an electron experiences is contaminated with 
self-interaction. Consequently, the negative of the energy eigenvalue of the highest occupied 
orbital in these approximate schemes often underestimates the first ionization potential. The 
HOMO-LUMO gap is also underestimated typically by about 30-40% 
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221 ] due to the so 



called band gap problem 23l, l24j|. In the HF theory, self-Coulomb potential exactly cancels 
the self-exchange potential for the occupied electrons and hence in the HF model the occupied 
electrons experience only electronic interactions with the other occupied electrons. Thus HF 
is self-interaction free for occupied orbitals. The virtual orbitals, however, experience the 
full electronic interaction with all the occupied electrons. Therefore, the HOMO-LUMO gap 
in HF theory is usually overestimated. 

The unscreened SOS method (Eq. Hj) is often used within the HF for qualitative purposes. 
A better and more appropriate expression in place of the HF eigenvalue differences in the 
denominator can be obtained based on the following argument. Once the zero-field ground- 
state HF determinant is found, it is possible to develop a multiconfigurational treatment 
which is successively more accurate by first constructing one-electron excitations, then two- 
electron excitations etc. Brillouin's theorem tells us that there is no coupling between 
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one-electron excitations and the occupied, zero-field molecular orbitals. When the many- 
electron field-induced perturbation F-ri) is added to this Hamiltonian, the only coupling 
between different zero-field single-electron excitations is due to this term and it is of the form 
(ip m \F ■ r\ip n ) for the excitation which replaces the formerly occupied state ip m by ip n . The 
energy of this state relative to the ground state is e m — e n — A mn with A mn a gap reducing 
correction given by: 

A mn = (lp m 1p m \l/r 12 \lpn^n) ~ (VW0nlV r 12|VW0n) (5) 

In the above expression the ipm are to be viewed as a product state of the spatial wavefunction 
and spinor. The second term is zero for spin-flip excitations. So, if a CI approach, with only 
single excitations, is employed, it is easy to show from perturbation theory that the correct 
expression for the polarizability within this approximation is not identical to the normal sum 
of states method because there is a pair-dependent correction to the energy denominators. 
In other words, rather than Eq. (0J, the asos should be 
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While the single-electron CI approach is itself an approximation, it is a more proper definition 
for the unscreened approximation and differs from the standard SOS method used in HF. 

As the HOMO-LUMO gap is the lowest of many excitation energies that contribute to 
asos, the trend in asos in two types of models (HF and DFT) should be, in general, opposite 
to that observed for the HOMO-LUMO gaps. While the general trends of overestimation 
and underestimation of asos^ the HF method and DF based models are to be expected, a 
detailed comparison on a large set of molecules is necessary to determine if this trend in the 
two types of models always holds, which partially motivates the present work. 

Our interest in the SOS method (Eq. H]) is primarily due to the fact that it is computation- 
ally inexpensive to apply, particularly for large, symmetric molecules for which the applica- 
tion of the finite-field method becomes very expensive. The SOS (using Eq. 
been employed in the past to compute the polarizability of carbon fullerenes 
In these works, the calculated asos were scaled using a scaling factor to account for the 
screening effects. As the point group symmetry of molecules can be efficiently employed 
in asoscalculations, Eq. (j4j) offers a possibility of obtaining estimates of polarizability of 



HP ) indeed has 
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large symmetric molecules that are beyond the reach of the more accurate, but symmetry- 
breaking, FF method. Our goal is to estimate the polarizabilities of carbon fullerenes, 
quantum dots, and nanocrystalline materials containing about 1000-2000 atoms by the SOS 
method [25J. The present work is a detailed investigation to understand the behavior of 
a>sos within different computational schemes. 

We have chosen a large set of 142 molecules belonging to the extended G2 set of molecules 
whose geometries are well known. To facilitate the comparison between different single- 
particle methods, we choose the same set of geometries and basis sets in our polarizability 
calculations in all models. This eliminates some sources of discrepancy in the comparison 
of polarizabilities. The geometries of 142 molecules were first optimized using the B3LYP 
hybrid DFT model using the 6-311G** orbital basis set. 26| B3LYP is one of the most pop- 
ular hybrid functionals. It contains a weighted mixture of the Becke88 exchange functional, 
the Lee- Yang-Parr gradient corrected correlation, Vosko-Wilk-Nusair correlation functional, 
the local-density exchange functional, and HF exchange [l|, ^j. The mixing coefficients are 
determined empirically by fitting to atomization energies. The HF and B3LYP calcula- 
tions were performed using the GaussianOsj^ at the Army Research Laboratory while the 
PBE calculations were performed using the NRLMOL suites of code. The NRLMOL is a 
massively parallelized suite of codes developed at the Naval Research Laboratory for DFT 



calculations for molecules 



281 ] . The molecular orbitals in the NRLMOL are expressed as a 



sum of linear combinations of Gaussian orbitals and the polarizability is computed using the 
FF method by numerically differentiating the calculated dipole moment at different field 
strengths/directions. The electric field vaule are chosen in the step of 0.005 a.u. These 
values provide reliable estimate of the polarizability. For a recent review on the details of 
polarizability calculations using the NRLMOL code, we refer the interested reader to Ref. 



29l . It is well known that for accurate computation of polarizability, well chosen Gaussian 



basis set augmented with diffused functions is often necessary. In this work, we have used 
the 6-31H-+G(d,p) type of basis for all polarizability calculations. This basis was selected 
on two criteria: (i) it is large enough to provide good accuracy for the polarizability (ii) its 
availability for all elements that belong to the chosen set of molecules. 

The asos an d ctpFof 142 molecules calculated within HF, B3LYP, and PBE models are 
compared in Table [H These are the mean polarizabilities obtained from the trace of the 
polarizability tensor as a = | Y^h=i a u- Although the main goal of this work is to compare 
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screened and unscreened polarizabilities, for completeness we also included in the table 
the experimental values of polarizabilities for selected molecules. All experimental values 
are from the CRC Handbook 30| . The overall agreement between theory and experiment is 
good. There are large differences for some molecules (e.g. C1F 3 which has a large vibrational 
contribution to polarizability). There are many possible causes for these differences such 
as temperature effects, basis sets effects and vibrational contributions. Proper treatment of 
these effects will give better agreement but is outside the scope and purpose of the present 
paper. Hence we shall not elaborate on the comparison of theoretical and experimental 
values of polarizability. The comparison of the (unscreened) asos with screened cuff, 
within the three models show different trends. In the HF method the asos are smaller than 
the o>ff with the exceptions of the P 2 , CS, COS, and N 2 molecules. The two values of 
polarizability for these exceptions are quite close, with the difference for P 2 being the largest. 
This almost always correct trend of HF polarizabilities is reversed in the parameter-free PBE 
DFT model. We find no exceptions. In PBE, SOS polarizabilities are always less than FF 
polarizabilities. The asos in PBE are consistently higher than the aFF polarizabilities. 
These trends, as noted earlier, are consistent with the general trend of the HF and PBE 
models to underestimate and overestimate the HOMO-LUMO energy gap. The asos are 
also predicted to be larger than aFF in the B3LYP method. However, the overestimation 
in the B3LYP model is less than that in the PBE model. The mean absolute error (MAE) 
calculated as MAE = J2l=i \ a sos~ ®ff\ is reduced from 3.55 A 3 for the PBE to 1.95 A 3 
for the B3LYP functional. The MAE is least in the HF theory, for which the mean error is 
negative. The hybrid B3LYP functional benefits from the opposite effects of overestimation 
in pure DF models and underestimation in the HF theory. 

TABLE I: The polarizabilities for the set of 142 (extended G2 set) of molecules 
computed by sum over states and coupled HF methods or finite field within 
the Hartree-Fock approximation, hybrid B3LYP model and the PBE gener- 
alized gradient approximation. The polarizabilities are calculated using the 
6-311++G(d,p) at the geometries optimized at B3LYP/6-311G** level. All 
values are in A 3 . I: SOS , II: Coupled HF or finite-field 

Molecule HF(I) HF(II) B3LYP(I) B3LYP (II) PBE (I) PBE (II) Expt.f30] 



H 2 


0.3 


0.4 


0.5 


0.4 


0.6 


0.5 


0. 


LiH 


1.9 


3.2 


4.0 


4.2 


5.5 


4.7 




BeH 


2.2 


4.2 


4.4 


4.5 


3.2 


4.7 




CH 


1.1 


1.7 


2.0 


1.8 


2.6 


1.8 
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CH 2 ( 3 Bi) 


1.1 


1.5 


1.8 


1.6 


1.6 


1.7 




CHaOAi) 


1.4 


1.9 


2.5 


2.1 


3.7 


2.1 




CH 3 


1.4 


1.8 


2.2 


2.0 


2.2 


2.0 




CH 4 


1.7 


2.1 


2.6 


2.2 


3.0 


2.2 


2.6 


NH 


0.7 


1.1 


1.2 


1.1 


1.5 


1.2 




NH 2 


1.0 


1.3 


1.6 


1.4 


2.0 


1.5 




NH a 


1.3 


1.6 


2.1 


1.8 


2.4 


1.8 


2.8 


OH 


0.5 


0.7 


0.8 


0.8 


1.0 


0.8 




H 2 


0.8 


1.0 


1.2 


1.1 


1.5 


1.1 




HF 


0.4 


0.5 


0.6 


0.5 


0.7 


0.6 




Li 2 


17.7 


30.3 


37.9 


29.3 


57.0 


30.2 


32.8 


LiF 


0.6 


0.8 


1.3 


1.2 


1.7 


1.4 




C 2 H 2 


2.6 


2.8 


4.3 


2.9 


5.1 


2.9 




C 2 H4 


3.1 


3.8 


5.1 


3.8 


6.1 


3.8 




C 2 Hg 


3.1 


3.8 


5.0 


4.0 


5.9 


4.1 




CN 


1.6 


1.9 


4.3 


2.9 


5.9 


2.7 




HCN 


2.0 


2.1 


3.3 


2.1 


3.9 


2.2 




CO 


1.4 


1.6 


2.5 


1.7 


3.0 


1.8 


1.9 


HCO 


1.7 


2.1 


3.0 


2.2 


3.7 


2.3 




H 2 CO (formaldehyde) 


1.8 


2.2 


3.1 


2.3 


3.8 


2.4 


2.8 


H3COH 


2.1 


2.6 


3.4 


2.8 


4.1 


2.9 




N 2 


1.4 


1.5 


2.4 


1.5 


2.8 


1.6 


1.7 


N 2 H 4 


2.3 


2.8 


3.8 


3.1 


4.5 


3.2 




o 2 


1.1 


1.4 


1.8 


1.2 


3.2 


1.2 


1.6 


H 2 2 


1.4 


1.8 


2.4 


1.9 


2.9 


1.9 




F 2 


0.7 


1.0 


1.3 


0.9 


1.6 


0.9 




co 2 


1.9 


2.0 


3.4 


2.2 


4.2 


2.2 


2.9 


SiH 2 ( 1 A 1 ) 


3.2 


4.4 


6.0 


4.5 


7.9 


4.6 




SiH 2 ( 3 Bi) 


2.6 


3.7 


4.7 


4.0 


4.3 


4.1 




SiH 3 


2.9 


3.9 


5.1 


4.2 


5.0 


4.3 




S1H4 


2.9 


3.7 


5.0 


4.0 


6.2 


4.2 


5.4 


PH 2 


2.5 


3.3 


4.4 


3.4 


5.3 


3.5 




PH 3 


2.9 


3.6 


5.0 


3.7 


6.0 


3.7 


4.8 


H 2 S 


2.2 


2.6 


3.6 


2.8 


4.3 


2.8 


4.0 


HC1 


1.3 


1.5 


2.0 


1.6 


2.3 


1.6 


2.8 


Na 2 


24.8 


40.2 


48.9 


34.0 


70.4 


35.4 




Si 2 


9.2 


13.8 


18.4 


14.6 


26.1 


19.1 




P 2 


7.1 


6.3 


13.1 


6.2 


16.7 


6.2 




s 2 


4.6 


5.5 


8.4 


4.8 


15.6 


4.7 




Cl 2 


2.8 


3.3 


4.7 


3.2 


5.7 


3.2 


4.6 


NaCl 


2.2 


2.9 


4.5 


4.0 


6.0 


4.4 




SiO 


3.2 


3.7 


6.0 


3.9 


7.7 


4.0 




cs 


3.5 


3.4 


6.4 


3.6 


8.1 


3.6 




so 


2.3 


2.8 


4.1 


2.7 


7.0 


2.7 




CIO 


1.7 


2.0 


3.5 


2.3 


4.5 


2.3 




C1F 


1.5 


1.8 


2.6 


1.9 


3.2 


1.9 
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Si 2 H 6 


6.3 


7.8 


CH3CI 


2.9 


3.4 


H 3 CSH 


3.9 


4.5 


HOC1 


2.0 


2.5 


S0 2 


3.1 


3.2 


BF3 


1.4 


1.7 


BCI3 


5.7 


6.2 


AIF3 


1.6 


2.0 


AICI3 


6.0 


6.9 


CF 4 


1.8 


2.1 


CCU 


7.6 


8.4 


COS 


4.1 


4.0 


cs 2 


7.4 


6.9 


COF 2 


1.8 


2.1 


S1F4 


1.9 


2.3 


Sicu 


7.9 


8.7 


N 2 o 


2.5 


2.4 


C2C14 


8.8 


10.0 


CF3CN 


3.5 


3.8 


CH3CCH (propync) 


4.1 


4.6 


CH2CCH2 (allcnc) 


4.6 


5.6 


C3H4 (cyclopropyne) 


4.3 


4.8 


CH3CHCH2 (propene) 


4.7 


5.5 


C 3 H 6 


4.4 


5.0 


C3H8 (propane) 


4.7 


5.5 


CH2CHCHCH2 (butidene) 


6.6 


8.0 


C4H6 (butyne) 


5.6 


6.5 


C4H6 (methylene cylcopropanc) 


5.9 


7.0 


C4H6 (bicyclobutane) 


5.7 


6.2 


C4H6 (cyclobutcne) 


5.8 


6.5 


C4H8 (Cyclobutane) 


5.8 


6.6 


C4H8 (isobutcnc) 


6.2 


7.2 


C4H10 (butane) 


6.2 


7.2 


C4H10 (isobutanc) 


6.2 


7.2 


CsHg (spiropentane) 


7.1 


8.1 


CeHg (benzene) 


9.4 


9.6 


CH2F2 (difiuromethylene) 


1.6 


2.0 


CHF3 (trifiuromcthane) 


1.7 


2.1 


CH2CI2 (dichloromethanc) 


4.3 


5.0 


CHCI3 (chloroform) 


5.9 


6.7 


CH3NH2 (methylamine) 


2.7 


3.2 


CH3CN (methyl cynaidc) 


3.4 


3.8 


CH3NO2 (nitromcthanc) 


3.7 


4.2 


CH3ONO (methyl nitrite) 


3.6 


4.0 


CHsSiHs (methyl silane) 


4.5 


5.5 


CHOOH (formic acid) 


2.3 


2.6 



11.4 


8.5 


14.2 


8.9 




4.6 


3.6 


5.4 


3.6 


5.3 


6.4 


4.8 


7.6 


4.8 




3.4 


2.5 


4.2 


2.5 




6.1 


3.4 


7.9 


3.5 


3.7 


2.4 


2.0 


3.0 


2.2 


3.3 


9.9 


7.0 


12.2 


7.2 


9.4 


3.0 


2.6 


3.8 


2.9 




10.7 


8.0 


13.1 


8.3 




3.0 


2.5 


3.6 


2.6 


3.8 


13.4 


9.1 


16.5 


9.4 


11.2 


7.7 


4.2 


9.7 


4.3 




14.7 


7.0 


19.2 


7.0 


8.7 


3.3 


2.4 


4.0 


2.5 




3.4 


2.9 


4.2 


3.1 


5.5 


13.9 


9.9 


17.0 


10.2 




4.8 


2.5 


6.1 


2.6 


3.0 


16.1 


10.8 


20.2 


11.1 




6.1 


4.3 


7.6 


4.5 




6.9 


4.9 


8.4 


5.0 


6.2 


8.0 


5.6 


9.9 


5.6 




7.0 


4.9 


8.4 


5.0 


6.2 


7.7 


5.7 


9.3 


5.8 




7.1 


5.3 


8.3 


5.4 




7.5 


5.8 


8.9 


6.0 


6.3 


11.9 


8.0 


15.1 


8.1 




9.6 


7.0 


11.7 


7.2 




10.0 


7.2 


12.1 


7.3 




9.3 


6.5 


11.0 


6.7 




9.6 


6.8 


11.6 


6.9 




9.3 


7.0 


11.0 


7.2 




10.4 


7.6 


12.5 


7.8 




10.0 


7.7 


11.8 


7.9 


8.2 


10.1 


7.7 


12.0 


7.9 


8.1 


11.7 


8.5 


13.9 


8.8 




15.6 


9.8 


18.9 


9.9 


10.0 


2.7 


2.3 


3.3 


2.4 




2.9 


2.4 


3.5 


2.5 


3.5 


7.2 


5.3 


8.7 


5.5 


6.5 


10.2 


7.2 


12.5 


7.5 


9.5 






o.r 


o.u 




5.7 


4.0 


6.9 


4.1 




6.9 


4.5 


8.7 


4.6 


7.4 


6.7 


4.5 


8.5 


4.7 




7.6 


5.9 


9.3 


6.2 




4.1 


3.0 


5.2 


3.2 


3.4 
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HCOOCH3 (methyl formate) 


3.7 


4.3 


6.4 


4.7 


7.8 


4.9 


5.0 


CH3CONH2 (acetamide) 


4.2 


4.9 


7.4 


5.5 


9.2 


5.8 


5.7 


C2H4NH (aziridine) 


3.9 


4.4 


6.3 


4.7 


7.5 


4.8 


2.6 


CNCN (cyanogen) 


4.1 


4.2 


8.0 


4.4 


10.4 


4.5 




(CH 3 ) 2 NH (dimethylaminc) 


4.2 


5.0 


6.9 


5.4 


8.3 


5.7 




CH3CH2NH2 (trans ethyalmine) 


4.2 


5.0 


7.0 


5.5 


8.4 


5.7 




CH2CO (ketene) 


3.2 


3.6 


5.8 


3.8 


7.1 


3.9 




C2H4O (oxirane) 


3.2 


3.7 


5.3 


4.0 


6.4 


4.1 




CH3CHO (acetaldehyde) 


3.3 


3.9 


5.7 


4.3 


7.0 


4.5 




HCOCOH (glyoxal) 


3.5 


4.0 


6.3 


4.4 


8.0 


4.6 




CH3CH2OH (cthanol) 


3.6 


4.3 


5.9 


4.7 


7.0 


4.8 




(CH 3 ) 2 (dimethylethcr) 


3.6 


4.3 


5.9 


4.7 


7.1 


4.9 




C2H4S (thioxirane) 


5.4 


5.9 


9.1 


6.1 


10.9 


6.2 




(CH 3 ) 2 SO (dimethyl sulfoxide) 


6.1 


6.9 


10.9 


7.6 


13.5 


7.9 




CH3CH2SH (cthanethiol) 


5.5 


6.4 


9.1 


6.7 


10.8 


6.9 




(CH 3 ) 2 S (dimethyl sulphide) 


5.5 


6.4 


9.2 


6.8 


11.0 


6.9 




CH2CHF (vinyl fluride) 


3.0 


3.6 


5.2 


3.8 


6.3 


3.9 




CH3CH2CI (ethyl chloride) 


4.5 


5.3 


7.4 


5.6 


8.8 


5.7 




CH2CHCI (vinyl chloride) 


4.4 


5.2 


7.6 


5.4 


9.3 


5.5 




CH3CHCN (acrylonitrile) 


5.1 


5.8 


9.0 


5.9 


11.4 


6.0 




(CH 3 ) 2 CO (acetone) 


4.8 


5.5 


8.1 


6.0 


9.9 


6.3 


6.3 


CH 3 COOH (acetic acid) 


3.7 


4.3 


6.4 


4.7 


7.8 


5.0 


5.1 


CH3COF (acetyl fluride) 


3.2 


3.7 


5.5 


4.1 


6.8 


4.3 




CH3COCI (acetyl chloride) 


4.9 


5.5 


8.7 


6.0 


10.8 


6.2 




CH3CH2CH2CI (propyl chloride) 


6.0 


7.0 


9.9 


7.5 


11.8 


7.7 




(CH 3 ) 2 CHOH (isopropanol) 


5.1 


6.0 


8.5 


6.5 


10.2 


6.8 




CH 3 CH 2 OCH 3 (methyl ethylether) 


5.0 


6.0 


8.3 


6.6 


10.0 


6.8 




(CH 3 ) 3 N (trimethylaminc) 


5.8 


6.7 


9.7 


7.4 


11.8 


7.8 




C4H4O (furan) 


6.2 


6.6 


10.5 


6.8 


12.7 


7.0 




C4H4S (thiophene) 


8.6 


8.7 


14.8 


9.0 


18.1 


9.1 


9.7 


C4H4NH (pyrole) 


7.0 


7.4 


11.6 


7.7 


14.0 


7.8 




C5H5N (pyridine) 


8.6 


8.7 


14.5 


9.1 


17.6 


9.2 


9.5 


CCH (ethynyl radical) 


2.7 


3.3 


4.9 


3.6 


5.9 


3.7 




CH 3 CO 


3.1 


3.8 


5.6 


4.2 


6.6 


4.5 




CH 2 OH (hydroxymcthyl) 


1.9 


2.4 


3.3 


2.8 


4.1 


2.9 




C1NO 


3.9 


4.4 


8.7 


4.3 


12.6 


4.4 




NF 3 


1.8 


2.2 


3.3 


2.4 


4.3 


2.5 


3.6 


PF 3 


2.5 


3.0 


4.6 


3.4 


5.8 


3.6 




o 3 


2.5 


2.7 


5.5 


2.4 


7.8 


2.4 


3.2 


F 2 


1.4 


1.8 


2.6 


1.8 


3.4 


1.8 




C1F 3 


3.0 


3.7 


6.4 


3.8 


8.8 


3.9 


6.3 


C2F4 


3.0 


3.5 


5.4 


3.9 


6.9 


4.1 




CH 3 (mcthoxy radical) 


1.9 


2.4 


3.3 


2.7 


4.2 


2.9 




CH 3 CH 2 


3.4 


4.1 


5.8 


4.6 


7.3 


4.8 




CH3S (mcthylsulfidc radical) 


3.3 


4.0 


5.5 


4.3 


6.7 


4.4 




CH3CH2 (ethyl radical) 


2.9 


3.5 


4.7 


3.9 


5.6 


4.0 
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As mentioned earlier, empirical schemes that convert the SOS polarizabilities to screened 



polarizabilities have been devised and applied to carbon fullerenes|13|, [14|, [16|, [31 1. These 
schemes were usually obtained in the random phase approximation and were applied to 
spherical molecules like fullerenes and cylindrical systems like nanotubes. Their application 
to molecules with arbitrary geometry has not yet been reported. An alternative possibility to 
improve upon the SOS polarizability will be to improve the approximation for the exchange- 
correlation functional. Correcting for the self-interaction in the approximate functional may 
improve the SOS polarizabilities. It has been found that the self-interaction corrections 
(SIC) significantly improves the eigenvalues and the band gapf^. The SIC implementation 
is however quite complicated to implement and computationally expensive [33]. The time 
dependent DFT (TD-DFT) within the linear response regime is another popular method 
for obtaining excitation energies but like the finite-field method it breaks the point-group 
symmetry. Much simpler approach would be to assume that the orbitals or the transition 
dipole moments obtained in the approximate DFT models to be reasonably accurate and 
correct only the eigenvalues using simple schemes. 

There are a number of methods in the literature that try to improve upon the HOMO- 
LUMO energy gap in the approximate DFT models typically by correcting the eigenvalues 
of the single-particle equations. Such corrections should result in some improvements in the 
SOS polarizability. Most of these methods are, however, quite complex and computationally 
expensive, limiting their advantages for correcting asos over the finite field method in 
terms of computational time. Hence simple schemes for correcting the eigenvalues are more 



appealing [34|, |35|, |36( and worth exploring in calculation of asos- Such simple schemes have 



been used in calculations of photoelectron spectra with some success [37]] . The simplest such a 
correction would to be replace the Ahl by the quasiparticle gap. The latter set of quantities 
can be computed by finite difference method or the so called ASCF method. This requires 
two additional self-consistent calculations, one for cation and other for anion. The correction 
5 to the A HL is then E(N - 1) + E(N + 1) - 2E(N) - e HOM o + e LUM o, where, E(N) is 
the self-consistent total energy of the system containing N electrons, enoMO/LUMO are the 
eigenvalues of the HOMO /LUMO of the neutral molecule. This correction, when applied to 
all single particle energies, would shift the occupied and unoccupied eigenvalues in opposite 
directions and hence may be useful in getting better estimate of asos- Its application to the 
L%2 for which the difference between the off and asos is large, indicate that this correction 
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overcorrects asos , roughly by a factor of two (from 57.0 A 3 to 19.0 A 3 for PBE). 

Another possible scheme can be devised by noting that the HF overestimates the asos 
whereas the PBE underestimates. It is possible to mix the HF asos an d the PBE asos to 
obtain finite-field quality polarizability. In absence of more accurate data, we choose B3LYP 
finite-field polarizabilities as a target set. The following interpolation scheme can be used: 
a mix = (1 — x) * affp + xaff^ s E , where the parameter x (0 < x < 1) can be determined 
by minimizing the mean absolute error (MAE) in the a mtx and a^ LY p- The minimization 
procedure gives the optimal mixing parameter to be 0.22. The MAE at the minimum is 0.4 
A 3 . Thus about 80% of HF asos mixed with about 20% of the PBE asos will give asos 
comparable to the B3LYP off- In this application more HF than DFT is required, whereas 
Becke needed less HF than DFT for optimizing B3LYP for atomization energies l] . This 
procedure requires calculation of asos i n two models which could be performed efficiently 
by making use of any symmetry that system may possess. Thus it could be applied to 
symmetric quantum dots, fullerenes and nanocrystals to obtain estimates of polarizability. 
We note that similar idea was explored for a few molecules by Dunlap and Karna[38]. These 
authors used average values of eigenvalues in the HF and the local-density approximation 
(DFT) eigenvalues in the SOS expression (Eq. Hj). 

In conclusion, a systematic comparison of the SOS and FF methods for the calculation 
of molecular dipole polarizability for a set of 142 molecules from the extended G2 set has 
been performed. The trends in the two sets, asos an d apF, of polarizabilities are examined 
in three widely used single-particle methods: The HF approximation, the PBE-GGA within 
DFT, and the hybrid B3LYP model that mixes DFT exchange with HF exchange. In order 
to minimize the other sources that can lead to differences in the polarizabilites in different 
models, the same set of molecular geometries (optimized at the B3LYP/6-311G** level) and 
6-311++G(d,p) orbital basis set was used. The calculations show that asos polarizabilities 
are almost always underestimated in the HF method. However, exceptions to the trend 
do exist. The asos of P2, CS, COS, and N 2 are overestimated with respect to the apF 
polarizability. On the other hand, in the PBE-GGA model, asos is always overestimated 
with respect to the off- These observations correlate with the generally observed trend of 
the respective underestimation and overestimation of the HOMO-LUMO gap in the PBE- 
GGA and HF models. Although the asos overestimates the apF m the hybrid B3LYP, 
the differences in the two polarizabilities is less than that observed for the PBE model. 
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The comparison of screened polarizabilities olff in the three models shows that the olff 
in B3LYP and PBE models are larger than those calculated with the HF method. Thus 
inclusion of correlation effects, in general, leads to increase in (finite-field) polarizability. 
Finally, a simple scheme that interpolates asos values in HF and PBE to obtain B3LYP 
ciff is devised by minimizing the mean absolute error. A simple scheme like this may be 
useful to estimate polarizabilities of large symmetric molecules such as fullerenes, quantum 
dots or nanocrystalline materials. 

Authors acknowledge Dr. S. P. Kama and Prof. P. B. Allen for discussions. The Gaus- 
sian03 calculations were performed at the Army Research Laboratory Major Shared Resource 
Center (ARL-MSRC). This work was supported in part by the Office of Naval Research, di- 
rectly and through the Naval Research Laboratory (ONR Grant No. N000140211046, and 
05PR07548-00) and in parts by the NSF through CREST grant (Grant No. NIRT-0304122), 
by the University of Texas at El Paso. 
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